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Abstract. - We measure the self-diffusion of colloidal rod-like virus fd in an isotropic and 
nematic phase. A low volume fraction of viruses are labelled with a fluorescent dye and dis- 
solved in a background of unlabelled rods. The trajectories of individual rods are visualized 
using fluorescence microscopy from which the diffusion constant is extracted. The diffusion 
parallel (-Dii) and perpendicular (D±) to the nematic director is measured. The ratio (D«/D±) 
increases monotonically with increasing virus concentration. Crossing the isotropic-nematic 
phase boundary results in increase of D« and decrease of D±_ when compared to the diffusion 
in the isotropic phase (Di so ). 



Introduction. Suspensions of semi-flexible polymers exhibit a variety of dynamical 
phenomena, of great importance to both physics and biology, that are still only partially 
understood. Advances over the past decade include direct visual evidence for a reptation-like 
diffusion of individual polymers in a highly entangled isotropic solution and shape anisotropy 
of a single polymer [1-4]. If the concentration of the polymers is increased, a suspension 
undergoes a first order phase transition to a nematic phase, which has long range orientational 
order but no long range positional order. As a result of the broken orientational symmetry 
it is expected that the diffusion of polymers in the nematic liquid crystals will be drastically 
different from that in concentrated isotropic solutions. While the static phase behavior of semi- 
flexible nematic polymers is well understood in terms of the Onsager theory and its extensions 
by Khoklov and Semenov [5,6], the dynamics of semi-flexible polymers in the nematic phase 
is much less explored [7] . 

In this paper, we determine the concentration dependence of the anisotropic diffusion of 
semi-flexible viruses in a nematic phase and compare it to the diffusion in the isotropic phase. 
Experimentally, the only data on the translational diffusion of colloidal rods in the nematic 
phase was taken in a mixture of labelled and unlabelled polydysperse boehmite rods using 
fluorescence recovery after photobleaching (FRAP) [8]. Theoretically, molecular dynamics 
simulations were performed on hard spherocylinder and ellipsoidal systems from which the 
anisotropic diffusion data was extracted [9-11]. The anisotropic diffusion has also been studied 
in low molecular weight thermotropic liquid crystals using NMR spectroscopy or inelastic 
scattering of neutrons [12]. 
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Real space microscopy is a powerful method that can reveal dynamics of colloidal and 
polymeric liquid systems that are inaccessible to other traditional techniques [3,7]. We use 
digital microscopy to directly visualize the dynamics of fluorescently labelled fd in a nematic 
background of the unlabelled fd. The advantage of this method is an easy interpretation of 
data and no need to obtain macroscopically aligned monodomains in magnetic fields. The 
advantages of using fd are its large contour length which can be easily visualized with optical 
microscope and its phase behavior which can be quantitatively described with the Onsager 
theory extended to account for electrostatic repulsive interactions and semi-flexibility [13,14]. 
Viruses such as fd and TMV have been used earlier to study the rod dynamics in the isotropic 
phase [15]. 

Experiment methods. - The physical characteristics of the bacteriophage fd are its length 
L=880 nm, diameter D=6.6 nm, persistence length of 2200 nm and a surface charge of 10 
e~/nm at pH 8.2 [16]. Bacteriophage fd suspension forms isotropic, cholcstcric and smec- 
tic phases with increasing concentration [16-18]. The free energy difference between the 
cholesteric and nematic phase is very small and locally the cholesteric phase is identical to ne- 
matic. We expect that at short time scales the diffusion of the rods for these two cases would 
be the same. Hereafter, we refer to the liquid crystalline phase at intermediate concentration 
as a nematic instead of a cholesteric. 

The fd virus was prepared according to a standard biological protocol using XLl-Blue 
strain of E. coli as the host bacteria [19]. The yields are approximately 50 mg of fd per liter 
of infected bacteria and virus is typically grown in 6 liter batches. Subsequently, the virus 
is purified by repetitive centrifugation (108,000 g for 5 hours) and re-dispersed in a 20 mM 
phosphate buffer at pH=7.5. First order isotropic-nematic (I-N) phase transition for fd under 
these conditions takes place at a rod concentration of 15.5 mg/ml. 

Fluorescently labelled fd viruses were prepared by mixing 1 mg of fd with 1 mg of suc- 
cinimidyl ester Alexa-488 (Molecular Probes) for 1 hour. The dye reacts with free amine 
groups on the virus surface to form irreversible covalent bonds. The reaction is carried out 
in small volume (100 /A, 100 mM phosphate buffer, pH=8.0) to ensure a high degree of la- 
belling. Excess dye was removed by repeated centrifugation steps. Absorbance spectroscopy 
indicates that there are approximately 300 dye molecules per each fd virus. Viruses labelled 
with fluorescein isothiocynante, a dye very similar to Alexa 488, exhibit the phase behavior 
identical to that of unlabelled virus. Since liquid crystalline phase behavior is a sensitive 
test of interaction potential, it is reasonable to assume that the interaction potential between 
labelled viruses is very similar to that between unlabeled viruses. 

The samples were prepared by mixing one unit of anti-oxygen solution (2 mg/ml glucose 
oxidase, 0.35 mg/ml catalase, 30 mg/ml glucose and 5% /3-mercaptoethanol) , one unit of a 
dilute dispersion of Alexa 488 labelled viruses and eight units of the concentrated fd virus 
suspension at the desired concentration. Under these conditions the fluorescently labelled 
viruses are relatively photostable and it is possible to continuously observe rods for 3-5 minutes 
without significant photoblcaching. The ratio of labelled to unlabelled particles is roughly kept 
at 1:30000. The samples were prepared by placing 4 fi\ of solution between a No 1.5 cover slip 
and coverslidc. The thickness of the samples is about 10 ^m. Thin samples are important to 
reduce the signal of out-of-focus particles. Samples are equilibrated for half an hour, allowing 
flows to subside and liquid crystalline defects to anneal. We have analyzed data at various 
distances from the wall and have not been able to observe a significant influence of wall on 
the diffusion of viruses. 

For imaging we used an inverted Nikon TE-2000 microscope equipped with 100 x 1.4 NA 
PlanApo oil immersion objective, a 100 W mercury lamp and a fluorescence cube for Alexa 
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Fig. 1 - (a) Image of fluorescently labelled rods dissolved in a background nematic phase of unlabelled 
rods. Scale bar is 5 /im. (b) Two-dimensional gaussian fit to a individual rod. Arrows indicate the 
long and short axis. The circle indicates the center of mass. From this fit it is possible to obtain the 
orientation of an individual fd rod. Pixel size is 129 nm. 

488 fluorescent dye. The images where taken with a cooled CCD Camera (CoolSnap HQ, 
Roper Scientific) set to an exposure time of 60 ms, running in a overlap mode at a rate of 16 
frames per second with 2x2 binning. The pixel size was 129 nm and the field of view was 89 
/ttm x 66 /im. Typically there were around hundred fluorescently labeled rods in the field of 
view. For each fd concentration ten sequences of four hundred images were recorded. 

Analysis method. - Figure shows a typical image of fluorescently labelled rods in a 
background nematic of unlabelled rods. Due to limited spatial and temporal resolution of the 
optical microscope, labelled fd appear as a slightly anisotropic rod, although the actual aspect 
ratio is larger then 100. To measure the anisotropic diffusion in the nematic phase, it is first 
necessary to determine the nematic director which has to be uniform within a field of view. 
Spatial distortion of the nematic would significantly affects our results. The centers of mass 
and orientation of rods are obtained sequentially. In a first step, a smoothed image is used 
to identify the rods and obtain the coordinates of its center of mass using image processing 
code written in IDL [20]. Subsequently, a two dimensional gaussian fit around a center of 
mass of each rod is performed (Fig-QJ)). From this fit the orientation of each rod-like virus is 
obtained. This procedure is then repeated for a sequence of images. 

The length of a trajectory is usually limited to a few seconds, after which the particles 
diffuse out of focus. In Fig. and b we plot the trajectories of an ensamble of particles 
for both isotropic and nematic sample. As expected the trajectories in the isotropic phase 
are spherically symmetric (Fig. while those in the nematic phase exhibit a pronounced 
anisotropy (Fig. |20- The symmetric nature of the distribution indicates that there is no 
drift or flow in our samples. We obtain the orientation of the nematic director using two 
independent methods. One method is to measure the main axis of the distribution shown 
in Fig. 12). This procedure assumes that the diffusion is largest along the nematic director. 
An alternative method is to plot a histogram of rod orientations which are obtained from 2D 
gaussian fits to each rod (Fig.^Jj). The resulting orientational distribution function (ODF) is 
shown in Fig. Et. In principle, it should be possible to obtain both the nematic director and 
order parameter from ODF shown Fig. El;. We find that the order parameter obtained in such 
a way is systematically higher then the order parameter obtained from more reliable x-ray 
experiment [14]. This is due to significant rotational diffusion each rod undergoes during an 
exposure time of 60 ms. 

The differences in the orientation of the nematic director obtained using these two methods 
is always less then 5 degrees. For the example shown in Figs. 03 and c, we obtain a nematic 
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Fig. 2 - (a) A collection of trajectories of fluorescently labelled virus particles in the isotropic phase. 
All trajectories are translated so that the first point is located at the origin. For clarity we only 
show the center of mass and not a line connecting subsequent point in a particle trajectory. The 
concentration of virus in this sample was 14 mg/ml. (b) Anisotropic trajectories of the fluorescently 
labelled viruses diffusing in the nematic phase. The concentration of the background virus in this 
sample was 21 mg/ml. x' and y' indicate a new lab-frame in which the director is aligned along the 
y' axis, (c) The orientational distribution function obtained by plotting the probability distribution 
function of the virus orientation for isotropic (open circles) and nematic phase (full squares). The 
orientation of the virus is obtained from two-dimensional gaussian fits, an example of which is shown 
in Fig. The angle of the nematic director obtained from (b) and (c) are almost identical. 



director at an angle of 31.2 ° while the peak of the orientational distribution function lies at 
30.2 °. The director can be "placed" along one of the two main axis by rotating the lab-frame. 

The diffusion coefficients of the rods parallel (D\\) and perpendicular (D±) to the director 
are calculated from the x'- and y '-component of the mean square displacement. When director 
lies along the y'-axis, D\\ and D± are given by: 

°ll = ]^£{24(*)-24(0)} 2 (i) 
£x = ^£{^(t)-^(0)} 2 , (2) 

where N is the number of traced particles. To obtain D j_, D x is multiplied with v2 since only 
one component of the diffusion perpendicular to the director is measured. The underlying 
assumption of our analysis is that the nematic director is oriented in the field of view. For 10 
fim thin samples this is reasonable. 

Results and discussion. - Typical mean square displacements (MSD) are shown in Fig. 
[3] for samples in an isotropic and nematic phase. On average the mean square displacement 
was linear over fifty frames in the nematic phase, but only about twenty five frames in the 
isotropic phase. The diffusion perpendicular to the director is slower in the nematic phase as 
compared to the isotropic phase. Therefore in the nematic phase, the particles stay longer in 
focus and can be tracked for a longer time. Since the MSD is linear over the entire time range 
and displacements are up to a few times the particle length, we are measuring pure long- 
time self-diffusion. Visual inspection of the trajectories in the concentrated isotropic phase, 
just below I-N coexistence shows no characteristics of the reptation observed in suspensions 
of long DNA fragments or actin filaments [2,3]. This points to the fact that fd is very 
weakly entangled in a concentrated isotropic suspension. This is in agreement with recent 
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Fig. 3 - The mean square displacement of rods along the director (full cubes) and perpendicular to the 
director (full triangles) for a nematic sample at virus concentration of 21 mg/ml. The isotropic data 
are given by the open points and were take just below I-N phase transition at virus concentration of 
14 mg/ml. The diffusion along the director is significantly enhanced when compared to the diffusion 
in the isotropic phase, while the diffusion perpendicular to the director is significantly suppressed. 
The mean square displacements shown in this figure are measured from a single field of view. 



micror neology measurements of fd suspensions [21]. We note that MSD's obtained from few 
hundred trajectories within a single field of view are very accurate. However, if we move to 
another region of the sample sample we obtain MSD with slightly different slope. This leads 
to conclusion that the largest source of error in measuring the anisotropic diffusion coefficient 
is the uniformity of the nematic director within the field of view. 

The concentration dependence of the anisotropic diffusion constants is shown in Fig. 0Ji. 
The nematic phase melts into a isotropic phase at low concentrations and freezez into a smectic 
phase at high concentrations. We made an attempt to measure the diffusion of rods in the 
smectic phase, but have not seen any appreciable diffusion on optical length scales over a 
time period of minutes. The most strinking feature of our data is a strong discontinuity in 
the behavior of the diffusion constant at the I-N phase transition. Compared to diffusions in 
isotropic case D iso , Dm is larger by a factor of four, while D± is smaller by a factor of two. 
The concentration dependence of Dm and D± exhibit different behavior. With increasing 
concentration, for Dm we measure an initial plateau, which is followed by a broad region 
where the diffusion rate decreases monotonically. Dx, however, shows a monotonic decrease 
of the diffusion constant over the whole concentration range where nematic phase is stable. 

It is useful to compare our results to previous theoretical and experimental work, especially 
the measurements of the diffusion coefficient for silica coated bochmitc rods [8] . In this work 
authors measure D»/D± w 2 for monodomain nematic samples which are in coexistence 
with isotropic phase. This is significantly different from D»/D± rs 7.5 for fd virus. Another 
significant difference is that results on bochmite indicate that both Dm and Dj_ are smaller 
then D lso . In contrast to our measurements where D\\ is much larger and D± is much smaller 
then D iso . 

When comparing our data to simulations of the diffusion of hard spherocylinders and el- 
lipsoids [9,11], one needs to compare equivalent samples. Scaling to rod concentration where 
the I-N transition takes place would be erroneous, since fd virus is a semi-flexible rod. The 
semi-flexibility of the virus drives the isotropic-nematic phase transition to higher concentra- 
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Fig. 4 - (a) The concentration dependence of the translational diffusion parallel to (-Dm) and per- 
pendicular to (D±) the nematic director are indicated by squares and triangles respectively. The 
nematic phase in coexistence with the isotropic phase occurs at c/d=15.5 mg/ml and is indicated by 
a vertical line. The x-axis is rescaled so that I-N transition takes place at [/d]jv = l. (b) The plot of 
the dimensionless ratio of the parallel to perpendicular diffusion constant D«/D perp as a function of 
the nematic order parameter. The concentration dependence of the nematic order parameter is taken 
from ref. [14]. Open triangles are data for hard spherocylinders with aspect ratio of 10 taken from 
ref. [11] while open circles are data for ellipsoids with aspect ratio 10 taken from [9] 



tions and it significantly decreases the order parameter of the nematic phase in coexistence 
with the isotropic phase [13,14]. We choose to compare data and simulations at the same 
value of the nematic order parameter which is determined independently [14]. For fd, the 
nematic order parameter is 0.65 at the I-N coexistence, monotonically increases with increas- 
ing rod concentration and saturates at high rod concentration. Experiment and simulation 
qualitatively agree and both show a rapid increase of Du/Dj_ ratio with increasing nematic 
order parameter (Fig.^J)). We note that there is a discrepancy between the simulations results 
obtain in references [9, 11] which might be due to different systems studied in these two paper. 

Interestingly, simulations predict that upon increasing rod concentration beyond I-N co- 
existence Du increases and subsequently upon approaching the smectic phase it will decrease. 
The authors argues that the non-monotonic behavior of is the result of the interplay 
between two effects. First, with increasing rod concentration the nematic order parameter 
increases which enhances D\\ . Second, with increasing rod concentration there is less free vol- 
ume which leads to decrease of Dm . The author further argues that the first effect dominates 
at low rod concentrations where the nematic order parameter rapidly increases while the sec- 
ond effect dominates at high rod concentrations where the nematic order parameter is almost 
saturated. In contrast, both of these effects contribute to a monotonic decrease in D± with 
increasing concentration, which is observed in simulations. Due to relatively large error in our 
experimental data, it is not clear if the behavior of Du is non-monotonic. There is an initial 
hesitation, but Du decreases over most of the concentration range. This difference between 
simulations and experiment might be because we compare experiments of semi-flexible fd to 
simulations of perfectly rigid rods. Compared to semi-flexible rods, the order parameter of 
rigid rods increases much faster with increasing rod concentration [14]. 

It would be of interest to extend our measurements to rotational diffusion in the isotropic 
and nematic phase. At present the rod undergoes significant rotational diffusion during each 
exposure which reduces resolution and prevents accurate determination of the instantaneous 
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orientation of a rod. It might be possible to significantly reduce the exposure time by either 
using a more sensitive CCD camera or a more intense laser as a illumination source. 

Conclusions. Using fluorescence microscopy we have visualized rod-like viruses and 
measured the anisotropic long-time self-diffusion coefficients in the isotropic and nematic 
phase. In the nematic phase the diffusion along the director and the diffusion perpendicular to 
the director decreases monotonically with increasing rod concentration. The ratio of parallel 
to perpendicular diffusion increases monotonically with increasing rod concentration. The 
results compare qualitatively with simulations on hard rods with moderate aspect ratios. 
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